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Abstract 



We study a stochastic multiplicative process with reset events. It is shown 
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(^ ■ that the model develops a stationary power-law probability distribution for 

On , the relevant variable, whose exponent depends on the model parameters. Two 

qualitatively different regimes are observed, corresponding to intermittent 

a I and regular behaviour. In the boundary between them, the mean value of 

the relevant variable is time- independent, and the exponent of the stationary 

n^ • distribution equals —2. The addition of diffusion to the system modifies in 

G . a non-trivial way the profile of the stationary distribution. Numerical and 

Q \ analytical results are presented. 
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The occurrence of power-law distributions (PLDs) is a common feature in the description 
of natural phenomena. These distributions appear in a wide class of nonequilibrium systems, 
ranging from physical processes such as dielectric breakdown, percolation, and rupture [|1|, 
to biological processes such as dendritic growth and large-scale evolution 0, to sociological 
phenomena such as urban development . Power-laws have been associated with the effect 
of the complex driving mechanisms inherent to these systems and with their intrincate 
dynamical structure. Criticality, fractals, and chaotic dynamics are known to be intimately 
related to PLDs 0]. 

In view of the ubiquity of PLDs in the mathematical description of Nature, much work 
has been recently devoted to detecting universal mechanisms able to give rise to such dis- 
tributions. In the frame of equilibrium processes, for instance, power-laws have been shown 
to derive from generalized maximum-entropy formulations ^ . For nonequilibrium phenom- 
ena, self-organized criticality (SOC) and stochastic multiplicative processes (SMPs) have 
been identified as sources of PLDs. According to the SOC conjecture [^], some nonequi- 
librium systems are continuously driven by their own internal dynamics to a critical state 
where, as for equilibrium phase transitions, power-laws are omnipresent. On the other hand, 
SMPs 1^ provide a (more flexible) mechanism for generating PLDs, based in the presence 
of underlying replication events. 

It is however well known that a pure SMP, 

n{t + l)=fi{t)n{t), (1) 

with n a random variable, does not generate a stationary PLD for n{t). Rather, it gives rise 
to a time-dependent log-normal distribution. To model the above mentioned phenomena, 
therefore, SMPs have to be combined with additional mechanisms. It has been shown that 



transport processes P , sources [Q , and boundary constraints |T0| are able to induce a SMP 
to generate power-laws. The aim of the present paper is to discuss an alternative additional 
mechanism, namely, randomly reseting of the relevant variable to a given reference value. 
In a real system, this would represent catastrophic annihilation or death events, seemingly 
originated outside the system. 

We consider a discrete-time stochastic multiplicative process n{t), added with reset events 
in the following way. At each time step, n is reset with probability g to a new value uq, drawn 
from a probability distribution PqIuq). If the reset event does not occur, n is multiplied by 
a random positive factor fi with probability distribution -P(yu). Namely, 

n(t + 1) = I ^^^^ ^ ^"^ ^^^'^ probability q, ,^. 

1 fi{t)n{t) with probability I — q. 

Between two consecutive reset events, n{t) thus behaves as a pure multiplicative process. 
When one of such events occurs, the multiplicative sequence starts again. 

In order to gain insight into the dynamics of process (0) we first consider the simplest 
case where no(t) and fi{t) are constant for all t. Since an arbitrary factor in the initial value 
of n is irrelevant to its subsequent evolution, we take uq = 1 without loss of generality. We 
have thus 

^C^ ^ 1) ^ I 1 with probability g, .^. 

1 fin{t) with probability I — q. 



This stochastic recursive equation can be readily solved to give 

, s _ j 1^^ with probability pk = q{l — q)^ {0 < k <t — 1), , . 

1 /i* with probability pt = (1 — g)*. 

Note that the possible values of n{t), Uk = fi^ {k = 0, l,...,t), lie in the interval [/x*, 1] 
for /i < 1 and in [1,/i*] for /i > 1. Except for the extreme value rit = /i*, the associated 
probabilities are time- independent. As time elapses, the probability of each possible value of 
n{t) is therefore quenched for n ^ fi'^, and the corresponding probability distribution evolves 
at this extreme value only. Thus, the distribution sequentially builds up in zones that lie 
increasingly further from n = 1. 

For large times, when the number of possible values of n{t) becomes also large, it is 
possible to give the probability distribution /(n) for n G (/i*, 1] for /i < 1 and n G [I,//*) for 
/i > 1 as a function of a continuous variable. In fact. 



f{n) = Pk 



dk 



dn 



ln/i| 



:n 



(5) 



with a = 1 — ln(l — q)/ In/i. In order to account for the contribution at ri = /i*, f{n) should 
be added with a delta-like term fo{t)5{n — /i*), where the factor /q can be obtained from the 
normalization of f{n). 

According to (|^), the stochastic process @ gives rise to a stationary ])ower- /aw distribu- 
tion f{n) in an increasingly large interval of values of n. For t -^ oo, f[n) is a stationary 
power-law distribution in (0, 1] for /i < 1, and in [1, oo) for ^ > 1. In contrast with multi- 
plicative processes with boundary constraints |T^ , there are no conditions on the parameters 
to obtain a stationary power-law distribution. For 1 — g < // < 1, the exponent of this distri- 
bution is positive (a < 0), and f{n) grows with n. In this situation, however, the distribution 
is defined for < ra < 1 and exhibits a cut off at n = 1. On the other hand, for fi < 1 — q 
or yU > 1 the exponent is negative [a > 0). For /i > 1, i.e. when n(t) G [1, oo), the moments 
rrii = J f{n)n'^dn diverge for z > a — 1, indicating the presence of intermittent amplifications 
||8|JTT|]. For n < 1 — q, rrii diverges ioi i < a — 1. 



It is interesting to relate the exponent of the power-law distribution with the evolution 
of the mean value {n{t)). From (^), this mean value can be written as 

<"W)^,3ii^.^l^f^.'(i-.)'. (a) 

For /i(l — g) < 1, the mean value of n(t) converges to a finite value (n) = g/[l — (1 — q)fj], 
whereas for /i(l — g) > 1 it "explodes". In the boundary between both regimes, where 
/i = 1/(1 — g), the exponent of the distribution is a = 2 and /(n) ~ n~^. This exponent is 
therefore to be associated with the explosion threshold. 

The power-law distribution in (|^) can also be inferred from a description of the evolution 
of fin). In fact, since at each time step where no reset occurs the probability contribution 
to f{n)dn comes from n' = n/fi, we can write 

ft+i{n)dn = (1 - q)ft{n/fi)d{n/iJ,). (7) 



Assuming now that this distribution is stationary, ft+i = ft, a. solution to (J^) is given by 
f{n) = An~°', with 

(1 - g)/i"-i = 1, (8) 

which produces the same value of a as in Eq. @. Note that (^ does not hold for n = 1, 
where the contributions to the probability come from reset events. 

The above argument provides a method for dealing with the general multiplicative process 
with reset events, Eq. (|^), when both fi and no are drawn from prescribed probability 
distributions P{fi) and Po^uq). We assume that Po('^o) is appreciably different from zero 
in a bounded region, where the contributions from reset events are relevant. Outside this 
region the evolution of f{n) can be written as 

POO 

ft+i{n)dn = {l-q) d^ P{n)ft{n/n)d{n/n), (9) 

Jo 

which generalizes Eq. ([^. Under the assumption of stationarity, this equation is solved by 
f{n) = An~", where the exponent a must verify 

(1-g)/ rf/x /i"-ip(/x) = 1. (10) 

For regular forms of PifJ^) this equation has at least one solution for a. When the probability 
is mainly concentrated in values of fi larger than unity the solution is expected to be positive 
(a > 0) and vice versa. 

As in the case of constant /i and tiq, a close relation exists here between the evolution of 
the average {n{t)) and the exponent of the power- law distribution. In particular, (n(t)) is 
found to remain stationary along the whole process when a = 2. Again, thus, the exponent 
a = 2 is associated with the explosion threshold, and marks the boundary between regular 
and intermittent evolution. This can be seen, for instance, from Eq. (^. Multiplication of 
this equation by n and integration over n yields 



(n(t + l)) = (l-g) 



dfi /iP(/i) 



{nit)). (11) 



Comparing with Eq. (|TUp, we readily note that the multiplicative constant (1 — (?) / dfi nP{ij) 
that governs the evolution of {n{t)) in Eq. (|lT]) equals unity for a = 2. 

In summary, depending on q and P{n) the system can be in a regular regime where 
{n{t)) converges to a finite value, or in an intermittence regime, where {n{t)) diverges. At 
the boundary, i.e. at the explosion threshold, {n{t)) remains constant and, independently of 
the specific value of q and of the particular form of -P(yu), the probability distribution f{n) 
exhibits a power-law tail with a characteristic exponent, f{n) ~ n~^. 

We have numerically checked that the exponent of the stationary profile of f{n) does not 
depend on the particular form of the distribution of reset values Po(^o)- Iii Fig- 1 we present 
the function f{n) obtained with constant /i and q, for three different choices of Po('^o)- A 
uniform distribution between no = and no = 1 (circles), an exponential distribution, 
-Po(^o) = (^o)~^ exp(— no/(no)), with (no) = 100 (squares), and a discrete distribution 
-Po(^o) = [^('^o ~ 1) +'^('^0 ■" 100)]/2. The particular form of Po sets a lower boundary for the 



region where f{n) behaves as a power law, but does not affect the corresponding exponent. 
Sohd hnes in the log- log plot of Fig. 1 have the theoretical slope a = 1.1054.... 

Figure 2 shows our simulation results for three different forms of P{fi)'- An exponential 
distribution P(/i) = (/i)~^exp(— /i/(/i)) with (fj,) = 2, a uniform distribution P(/i) = 5/2 
with /i G [9/10, 13/10], and a discrete distribution P(yu) = J2k=i ^(a* ~ yUfc)/3 with /ii = 1, 
fi2 = 6/5 and fi^ = 7/5. The slope of the solid lines has been obtained numerically for 
various values of q from Eq. (|TD|). This yields a = 1.4965... for the exponential distribution 
with q = 0.2, a = 1.2195... for the uniform distribution with q = 0.02, and a = 1.8965... for 
the discrete distribution with q = 0.15. In all cases, our numerical and analytical results are 
in full agreement within six to nine decades in the power-law region. 

We have also investigated the effects of diffusive transport on the process (|). With this 
aim, we have considered a one-dimensional array of elements whose individual dynamics is 
given by (|]) and, at each time step, we have incorporated an interaction mechanism that 
mimics diffusion. After the multiplicative process with reset events has been applied, the 
state of each element is further changed to 

n'M = (1 - D)n,{t) + - K+i(t) + n,_i(t)] , (12) 

where i labels the elements in the array, with periodic boundary conditions. Then, n[ is used 
as the input state for the next step. In this deterministic, time-discrete version of diffusive 
transport, D plays the role of a diffusion constant. 

Figure 3 summarizes our numerical results on the effect of diffusion on the SMP @, 
displaying the dependence of the power-law exponent with the diffusion constant. We have 
chosen values of q and fi such that the different regimes of the process have been explored. 
The value of the multiplicative constant has been fixed in this case to /z = 4/3. In the 
regular regime (i.e. /i(l — q) < 1), diffusion produces a decrease of a in the power-law 
distribution. This can be understood if we consider that the role of diffusion is to deplete 
dense areas, transporting material to less occupied cells. The multiplicative process is not 
fast enough in this regime to balance the joint effect of reset events and diffusion. As a 
result, underpopulation occurs in the high-density region, and a decreases (g = 0.3 in Fig. 
3). In the intermittent regime (g = 0.23 i.e. /i(l — q) > 1), diffusion favors the opposite 
effect. Remarkably, diffusion does not have any effect on the value of a when the system is 
evolving at the explosion threshold. Within numerical errors, in fact, a = 2 irrespectively 
of the value oi D. It is also worth to point out that the qualitative behaviour of the process 
depends on fj, and q only. Changing D does not allow the system to switch between the 
intermittent and the regular regimes. 

Summing up, in this paper we have studied a stochastic multiplicative process with reset 
events. The combination of this random reseting with the replication events driven by the 
stochastic process allows for the development of a stationary distribution in the system, both 
when the mean value of the relevant variable converges to a finite value (regular regime) 
and when it diverges (intermittent regime). The regime at the boundary between regular 
and intermittent behaviour is of particular interest. At this point, where the overall effects 
of the multiplicative process are exactly balanced by the random resets, the mean value of 
the relevant variable remains constant in time. We have shown that this property is closely 
related with the fact that the exponent of the power-law stationary distribution equals —2. 



This value is to be related with Zipf law, which predicts the same exponent of power-law 
distributions in a series of seemingly disparate natural systems |@,|^. Thus, the SMP with 
reset events offers an alternative explanation of this ubiquitous exponent. In fact, whereas 
a general trend of biological and social systems could be to improve their growth rates by 
increasing the parameter /i, it is on the other hand to be expected that external constrains 
are going to operate in order to avoid divergencies by increasing q. It is not unlikely that the 
competition between these two processes could lead real systems to this boundary between 
regular behavior and developed intermittency. 
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FIGURES 

FIG. 1. Stationary distribution f{n) for /i = 1.1 and q = 0.01, and for different distributions 
of reset values PoC^-o) (see text). Straight lines have the theoretical slope a = 1.2195 .... 

FIG. 2. Stationary distributions f{n) for different forms of i-*(/u) and different values of q (see 
text). The slope of the straight lines has been obtained through numerical solution of (p^. 

FIG. 3. Dependence of the exponent a on the diffusion coefficient D for // = 4/3 and three 
values of q corresponding to the intermittence regime {q = 0.23), the regular phase (g = 0.3), and 
the explosion threshold {q = 0.25). The error bars stand for the error of a in a least square fit to 
the numerical data. 
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